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Abstract 

We discuss in details the role of Wigner 6 j symbol as the basic building block uni- 
fying such different fields as state sum models for quantum geometry, topological 
quantum field theory, statistical lattice models and quantum computing. The appar- 
ent twofold nature of the 6j symbol displayed in quantum field theory and quantum 
computing -a quantum tetrahedron and a computational gate- is shown to merge 
together in a unified quantum-computational SU (2)-state sum framework. 

Keywords: quantum theory of angular momentum; Wigner Qj symbol; discretized 
quantum gravity; spin network quantum simulator 
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1 Introduction 




The above illustration shows a variant woodcut printer's device on verso last leaf of 
a rare XVI century edition of Plato's Timaeus, (Divini Platonis Operum a Marsilio 
Ficino tralatorum, Tomus Quartus. Lugduni, apud Joan Tornaesium M.D.XXXXX). 
The printer's device to the colophon shows a medaillon with a tetrahedron in cen- 
tre, and the motto round the border: Nescit Labi Virtus, Virtue cannot faifl. This 
woodcut beautifully illustrates the role of the perfect shape of the tetrahedron in 
classical culture. The tetrahedron conveys such an impression of strong stability as 
to be considered as an epithome of virtue, unfailingly capturing us with the depth 
and elegance of its shape. However, as comfortable as it may seem, this time- 
honored geometrical shape smuggles energy into some of the more conservative 
aspects of Mathematics, Physics and Chemistry, since it is perceptive of where the 
truth hides away from us: the quantum world. As Enzo says, the geometry of the 
tetrahedron actually takes us on a trip pointing to unexpected connections between 
the classical and the quantum. He has indeed often entertained us with descriptions 
of open terrains of Physics and Chemistry which are bumpy, filled with chemical 
bonds and polyhedra, and which bend abruptly in unexpected directions. We do 
feel that, like any good adventure, it is not the destination, but what we unexpect- 
edly found around the bend that counts. Thus, the story we wish to tell here is the 
story of what, together with Enzo, we found around the bend: the unfailing virtues 
of the quantum tetrahedron. 

'a more pedantic rendering is: Virtue ignores the possibility of sliding down. 
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Our story starts by recalling that the (re)coupling theory of many SU (2) angu- 
lar momenta -framed mathematically in the structure of the Racah-Wigner tensor 
algebra- is the most exhaustive formalism in dealing with interacting many-angular 
momenta quantum systems As such it has been over the years a common 

tool in advanced applications in atomic and molecular physics, nuclear physics as 
well as in mathematical physics. Suffices here to mention in physical chemistry 
the basic work of Wigner, Racah, Fano and others (see the collection of reprints 
and the Racah memorial volume quoted in Q below) as well as the recent book 
Hi on topics covered in this special issue. 

In the last three decades there has been also a deep interest in applying (exten- 
sions of) such notions and techniques in the branch of theoretical physics known 
as Topological Quantum Field Theory, as well as in related discretized models for 
3-dimensional quantum gravity. More recently the same techniques have been 
employed for establishing a new framework for quantum computing, the so-called 
"spin network" quantum simulator. 

In previous work in collaboration with Enzo (21 we have stressed the combina- 
torial properties of Wigner 6j symbols (and of its generalizations, the 3nj symbols, 
see O) which stand at the basis of so many different fields of research. 

The aim of the present paper is to discuss in details the apparent twofold na- 
ture of the 6j symbol displayed in quantum field theory and quantum computing, 
and to convey the idea that these two pictures actually merge together. In section 
2 the 6j is looked at as a real "tetrahedron", the basic magic brick in constructing 
3-dimensional quantum geometries of the Regge type, while in section 3 it plays 
the role of a magic box, namely the elementary universal computational gate in a 
quantum circuit model. Thus the underlying physical models embody, at least in 
principle, the hardware of quantum computing machines, while a quantum com- 
puter of this sort, looked at as a universal, multi-purpose machine, might be able 
to simulate "efficiently" any other discrete quantum system. More remarks this 
topic are postponed to the end of section 3, while most mathematical definitions 
and results on Wigner 6j symbols needed in the previous sections are collected in 
Appendix A. 

2 Tetrahedra and 6j symbols in quantum gravity 

From a historical viewpoint the Ponzano-Regge asymptotic formula for the 6j 
symbol [7], reproduced in (fl4l) of Appendix A.l, together with the seminal paper 
||8l in which "Regge Calculus" was founded, are no doubt at the basis of all "dis- 
cretized" approaches to General Relativity, both at the classical and at the quantum 
level. 

In Regge's approach the edge lengths of a "triangulated" spacetime are taken 
as discrete counterparts of the metric, a tensorial quantity which encodes the dy- 
namical degrees of freedom of the gravitational field and appears in the classical 
Einstein-Hilbert action for General Relativity through its second derivatives com- 
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bined in the Riemann scalar curvature. Technically speaking, a Regge spacetime 
is a piecewise linear (PL) "manifold" of dimension D dissected into simplices, 
namely triangles in D = 2, tetrahedra in D = 3, 4-simplices in D = 4 and so on. 
Inside each simplex either an Euclidean or a Minkowskian metric can be assigned: 
accordingly, PL manifolds obtained by gluing together D-dimensional simplices 
acquire an overall PL metric of Riemannian or Lorentzian signature^. 
Consider a particular triangulation T D (£) — > M D , where A4 D is a closed, locally 
Euclidean manifold of fixed topology and I denotes collectively the (finite) set of 
edge lengths of the simplices in T D '. The Regge action is given explicitly by (units 
are chosen such that the Newton constant G is equal to 1) 

S(T D {£)) = S D {£) = Vol (D - 2) (a,) e t , (1) 

where the sum is over (D — 2)-dimensional simplices cij 6 T D (called hinges 
or "bones"), Vol^ D ~ 2 ^ (<Tj) are their (D — 2)-dimensional volumes expressed in 
terms of the edge lengths and e t represent the deficit angles at Cj. The latter are 
defined, for each i, as 2ir — ^ fc 0^, where 6*^ are the dihedral angles between 
pairs of (D — l)-simplices meeting at cxj and labeled by some k. Thus a positive 
[negative or null] value of the deficit angle ej corresponds to a positive [negative or 
null] curvature to be assigned to the bone i, detected for instance by moving a D- 
vector along a closed path around the bone i and measuring the angle of rotation. 
Even such a sketchy description of Regge geometry should make it clear that a dis- 
cretized spacetime is flat (zero curvature) inside each D-simplex, while curvature 
is concentrated at the bones which represent "singular" subspaces. It can be proven 
that the limit of the Regge action (Q]) when the edge lengths become smaller and 
smaller gives the usual Einstein-Hilbert action for a spacetime which is "smooth" 
everywhere, the curvature being distributed "continuously". Regge equations -the 
discretized analog of Einstein field equations- can be derived from the classical 
action by varying it with respect to the dynamical variables, i.e. the set {£} of edge 
lengths of T D (£), according to Hamilton principle of classical field theory (we re- 
fer to [9] for a bibliography and brief review on Regge Calculus from its beginning 
up to the 1990's). 

Regge Calculus gave rise in the early 1980's to a novel approach to quantization 
of General Relativity known as Simplicial Quantum Gravity (see (9j [TUJ [TU and 
references therein). The quantization procedure most commonly adopted is the 
Euclidean path-sum approach, namely a discretized version of Feynman's path- 

2 Einstein's General Relativity corresponds to the physically significant case of a 4-dimensional 
spacetime endowed with a smooth Lorentzian metric. However, models formulated in "non- 
physical" dimensions such as D = 2, 3 turn out to be highly non trivial and very useful in a variety 
of applications, ranging from conformal field theories and associated statistical models in D = 2 to 
the study of geometric topology of 3-manifolds. Moreover, the most commonly used quantization 
procedure of such theories has a chance of being well-defined only when the underlying geometry is 
(locally) Euclidean, see further remarks below. 
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integral describing D-dimensional Regge geometries undergoing "quantum fluc- 
tuations" (in Wheeler's words a "sum over histories" Ifl2l . formalized for gravity in 
the so-called Hawking-Hartle prescription lfT3l ). Without entering into technical 
details, the discretized path-sum approach turns out to be very useful in addressing 
a number of conceptual open questions in the approach relying on the geome- 
try of smooth spacetimes, although the most significant improvements have been 
achieved for the D = 3 case, which we are going to address in some details in the 
rest of this section. 

Coming to the interpretation of Ponzano-Regge asymptotic formula for the 
6j symbol given in (PT4l ) of Appendix A. 1 , we realize that it represents the semi- 
classical functional, namely the semiclassical limit of a path-sum over all quan- 
tum fluctuations, to be associated with the simplest 3-dimensional "spacetime", 
an Euclidean tetrahedron T. In fact the argument in the exponential reproduces 
the Regge action S 3 (€) for T since in the present case (D — 2) simplices are 1- 
dimensional (edges) and Vo\( D ~ 2 \ai) in (Q} are looked at as the associated edge 
lengths, see the introductory part of Appendix A. 

More in general, we denote by T 3 (J) — > M 3 a particular triangulation of a 
closed 3-dimensional Regge manifold M 3 (of fixed topology) obtained by assign- 
ing SU(2) spin variables {j} to the edges of T 3 ■ The assignment must satisfy a 
number of conditions, better illustrated if we introduce the state functional associ- 
ated with T 3 (j), namely 



Z[T 3 U)^M 3 ;L]=A(Ly N °T](-l) 2 ^ WA y[ ( f )B P 1 ^ *H ( 2 ) 



where Nq, N\, N% are the number of vertices, edges and tetrahedra in T 3 (j), 
A(L) = 4L 3 /3C (L is a fixed length and C an arbitrary constant), = (2Ja + 1) 
are the dimensions of irreducible representations of SU (2) which weigh the edges, 

4>b = {— l) Zj p= i: ' p and {:::}b are 6j symbols to be associated with the tetrahedra 
of the triangulation. Finally, the Ponzano-Regge state sum is obtained by summing 
over triangulations corresponding to all assignments of spin variables {j} bounded 
by the cut-off L 



where the cut-off is formally removed by taking the limit in front of the sum. 

It is not easy to review in short the huge number of implications and further 
improvements of Ponzano-Regge state sum functional Q, as well as its deep and 
somehow surprising relationships with so many different issues in modern theoret- 
ical physics and in pure mathematics. We are going to present in the rest of this 
section a limited number of items, whose selection is made mainly on the basis 



A=l 



B=l 





(3) 



{j}<L 
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of their relevance for (quantum) computational problems raised in the next section 
(we remind however the importance of this model in the so-called "loop" approach 
to quantum gravity fl4l . see also iTTTlD . 

(a) As already noted in [7], the state sum Zpr [A4 3 ] is a topological invariant 
of the manifold A4 3 , owing to the fact that its value is actually indepen- 
dent of the particular triangulation, namely does not change under suitable 
combinatorial transformations. Remarkably, these "moves" are expressed 
algebraically in terms of the relations given in Appendix A.2, namely the 
Biedenharn-Elliott identity (PT71) -representing the moves (2 tetrahedra) f-> 
(3 tetrahedra)- and of both the Biedenharn-Elliott identity and the orthogo- 
nality conditions (fT8l) for Qj symbols, which represent the bary centric move 
together its inverse, namely (1 tetrahedra) <-> (4 tetrahedra). 

(b) In lfl5l a "regularized" version of (O -based on representation theory of a 
quantum deformation of the group SU{2)- was proposed and shown to be a 
well-defined quantum invariant for closed 3-manifold^. 

Its expression reads 



where the summation is over all {j} labeling highest weight irreducible rep- 
resentations of SU(2) q (q = exp{27ri/r}, with {j = 0, 1/2, 1 . . . ,r — 1}), 
wa = ( — 1 ) 2j A [2jA + 1] q where [ ] q denote a quantum integer, w = 2r/(q — 
q^ 1 ) 2 and | ::: \b represents here the q-6j symbol whose entries are the 
angular momenta ji,i = 1, . . . , 6 associated with tetrahedron B. If the de- 
formation parameter q is set to 1 one gets Ztv [-A4 3 ; 1] = T^pr [A4 3 ]. 
It is worth noting that the q-Racah polynomial -associated with the q-6j 
by a procedure that matches with what can be done in the SU (2) case, see 
([TBI in Appendix A.2- stands at the top of Askey's q-hierarchy collecting 
orthogonal q-polynomials of one discrete or continuous variable. On the 
other hand, the discovery of the Turaev-Viro invariant has provided major 
developments in the branch of mathematics known as geometric topology 



(c) The Turaev-Viro or Ponzano-Regge state sums as defined above can be gen- 
eralized in many directions. For instance, they can be extended to simpli- 
cial 3-manifold endowed with a 2-dimensional boundary ifTTl and to D- 
manifolds [18] (giving rise to topological invariants related to suitable (dis- 
cretized) topological quantum field theory of the Schwarz type |fl9l ). 

3 The adjective "quantum" refers here to "deformations" of semi-simple Lie groups introduced by 
the Russian School of theoretical physics in the 1980's in connection with inverse scattering theory. 
From the mathematical viewpoint the Turaev-Viro invariant, unlike the Ponzano-Regge state sum 
functional, is always finite and has been evaluated explicitly for some classes of 3-manifolds. 




(4) 



[16]. 
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(d) The fact that the Turaev-Viro state sum is a topological invariant of the un- 
derlying (closed) 3-manifold reflects a crucial physical property of gravity 
in dimension 3 which makes it different from the corresponding D = 4 case. 
Loosely speaking, the gravitational field does not possess local degrees of 
freedom in D = 3, and thus any "quantized" functional can depend only on 
global features of the manifold encoded into its overall topology Actually 
the invariant (0]) can be shown to be equal to the square of the modulus of the 
Witten-Reshetikhin-Turaev invariant, which in turn represents a quantum 
path-integral of an SU (2) Chern-Simons topological field theory -whose 
classical action can be shown to be equivalent to Einstein-Hilbert action 
[20]- written for a closed oriented manifold M 3 ETT [22Tl . Then there exists 
a corresponence 

Z TV [M 3 ;q] <— >■ \Z WRT [M 3 ;k]\ 2 , (5) 

where the "level" k of the Chern-Simons functional is related to the defor- 
mation parameter q of the quantum group. 

Despite the "topological" nature of Turaev-Viro (Ponzano-Regge) state sum and 
Witten-Reshetikhin-Turaev functionals in case of closed 3-manifolds, whenever 
a 2D-dimensional boundary occurs in M 3 , giving rise to a pair (Ai 3 , £), where 
£ is an oriented surface (or possibly the disjoint union of a finite number of sur- 
faces), things change radically. For instance, if we add a boundary to the man- 
ifold in Witten-Reshetikhin-Turaev quantum functional, the theory induced on 
£ is a Wess-Zumino-Witten (WZW)-type Conformal Field Theory (CFT) EOT 
endowed with non-trivial quantum degrees of freedom. In particular, the frame- 
works outlined above can be exploit to establish a direct correspondence between 
2D Regge triangulations and punctured Riemann surfaces, thus providing a novel 
characterization of the WZW model on triangulated surfaces on any genus 1231 at 
a fixed level k. 

We cannot enter here into many technical details on these developments. It should 
be sufficient to remark that, when addressing "boundary" CFT, the geometric role 
of the quantum tetrahedron shades out, while its algebraic content is enhanced 
given that the (q)-6j-symbol plays the role of a "duality" (or "fusion") matrix, 
similar to a "recoupling coefficient" between different basis sets, as (TTTT) in Ap- 
pendix A suggests. 

(e) In 11241 a (2 + l)-dimensional decomposition of Euclidean gravity (which 
takes into account the correspondence ((51)) is shown to be equivalent, under 
mild topological assumptions, to a Gaussian 2D fermionic system, whose 
partition function takes into account the underlying 3D topology. More pre- 
cisely, the partition function for free fermions propagating along "knotted 
loops" inside a 3-dimensional sphere corresponds to a 3D Ising model on 
so-called knot-graph lattices. On the other hand, the formal expression of 
3D Ising partition function for a dimer covering of the underlying graph 
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lattice can be shown to coincide with the permanent of the generalized in- 
cidence matrix of the lattice [25] [26). Recall first that the permanent of an 
n x n matrix A is given by 

n 

per[A] = IJ 'MO (6) 

where % j(T m are minors of the matrix, a(i) is a permutation of the index 
i = 1, 2, . . . , n and S n is the symmetric group on n elements. A graph lat- 
tice (S associated with a fixed orientable surfaces S of genus g embedded in 
S 3 may be constructed by resorting to the so-called "surgery link" presen- 
tation. Then the incidence matrix of such piecewise linear graph with, say, 
n vertices, is defined as an n x n matrix A = (a^) with entries in (1,0) 
according to whether vertices i, j are connected by an edge or not. Finally, 
the Ising partition function turns out to be a weighted sum -over all possible 
configurations of knot-graph lattices- of suitable "determinants" of general- 
ized forms of the incidence matrices which take into account the topology of 
the underlying manifold. We skip however other technical details and refer 
to G71 for a short account of these results (which will be briefly reconsidered 
in the following section in the context of quantum computational questions). 

The deep relationship between 3D quantum field theories that share a "topo- 
logical" nature and (solvable) lattice models in 2D, sketched in the last item by 
resorting to a specific example, was indeed predicted in the pioneering paper by 
E. Witten |[28l . Not so surprisingly, the basic quantum functional that realizes 
this connection was identified there with the expectation value of a certain tetra- 
hedral configuration of braided Wilson lines, where "Wilson lines" are quantum 
observables associated with "particle trajectories" that in general look like sheafs 
of braided strands propagating from a surface Si to another £2, both embedded in 
a 3D background. 

3 6j symbol and quantum algorithms 

The model for universal quantum computation proposed in ||29l , the "spin net- 
work" simulator, is based on the (re)coupling theory of SU(2) angular momenta 
as formulated in the basic texts HJ|2]| on the quantum theory of angular momentum 
and the Racah-Wigner algebra respectively. At the first glance the spin network 
simulator can be thought of as a non-Boolean generalization of the Boolean quan- 
tum circuit model^ lf30l . with finite-dimensional, binary coupled computational 

4 Recall that this scheme is the quantum version of the classical Boolean circuit in which strings 
of the basic binary alphabet (0, 1) are replaced by collections of "qubits", namely quantum states in 
(C 2 )® N , and the gates are unitary transformations that can be expressed, similarly to what happens 
in the classical case, as suitable sequences of "elementary" gates associated with the Boolean logic 
operations and, or, not. 
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Hilbert spaces associated with N mutually commuting angular momentum opera- 
tors and unitary gates expressed in terms of: 

i) recoupling coefficients (3nj symbols) between inequivalent binary coupling schemes 
of N = (n + 1) 5C/(2)-angular momentum variables (j'-gates); 

ii) Wigner rotations in the eigenspace of the total angular momentum J (M-gates) 
(that however will not be taken into account in what follows, see section 3.2 of 11291 
for details) 



BOOLEAN Q-CIRCUIT 



TOPOLOGICAL QFT 



Q-AUTOMATA 



GENERALIZED 
Q-CIRCUIT 



SPIN NETWORK 
Q-SfMULATOR 



STATE SUM 
MODELS 



In the diagram we try to summarize various aspects of the spin network sim- 
ulator together with its relationships with other models for Q-computation, in the 
light of underlying physical frameworks discussed in the previous section. 
On the left-hand portion of the diagram the standard Boolean quantum circuit 
is connected with a double arrow to the so-called topological approach to quan- 
tum computing developed in OTl (based, by the way, on the Witten-Reshetikhin- 
Turaev approach quoted in item (d) of the previous section). This means in prac- 
tice that these two models of computation can be efficiently "converted" one into 
the other. The Boolean case is connected one-way to the box of the generalized 
Q-circuit because it is actually a particular case of the latter when all N angular 
momenta are ^-spins. 

On the right-hand column, the double arrows stemming from the box of the spin 



9 



network Q-simulator relate it to its reference models: from the viewpoint of quan- 
tum information theory it is a generalized Q-circuit, as already noted before, while 
its physical setting can be assimilated to state sum-type models discussed in the 
first part of the previous section. 

The upper arrow is to be meant as generating, from the general Q-computational 
scheme, families of "finite-states" Q-automata able to process in an efficient way 
a number of specific algorithmic problems that on a classical computer would re- 
quire an exponential amount of resources (cfr. the end of this section). 

Besides the features described above, the kinematical structure of the Q-spin 
network complies with all the requisites of an universal Q-simulator as defined by 
Feynman in P2l . namely 

• locality, reflected in the binary bracketing structure of the computational Hilbert 
spaces, which bears on the existence of poly-local, two-body interactions; 

• discreteness of the computational space, reflected in the combinatorial structure 
of the (re)coupling theory of SU(2) angular momenta El l33l l34l : 

• discreteness of time, given by the possibility of selecting controlled, step-by-step 
applications of sequences of unitary operations for the generation of (any) process 
of computation; 

• universality, guaranteed by the property that any unitary transformation operat- 
ing on binary coupled Hilbert spaces (given by SU(2) 3nj symbols) can be recon- 
structed by taking a finite sequence of Racah-Wigner transforms implemented by 
expression of the type given in (TTTT) of Appendix A (possibly apart from phases 
factors), as shown in [2], topic 12. 

Then the Wigner 6j symbol plays a prominent role also in the spin network Q- 
simulator scheme, where it is the "elementary" unitary operation, from which any 
"algorithmic" procedure can be built up. The meaning of the identities (flTT ) (fT8T ) 
satisfied by the 6j's in the present context is analyzed at length in ||29l , (section 
4.2 and Appendix A) and can be related to the notion of intrinsic "parallelism" of 
quantum computers. 

A caveat is however in order: the complexity class of any classical [quantum] 
algorithm is defined with respect to a "standard" classical [quantum] model of com- 
putational At the quantum level, such a reference model is the Boolean Q-circuit 
ll30l . and thus what is necessary to verify is that a 6j symbol with generic entries 
can be efficiently (polynomially) processed by a suitably designed Q-circuit. Note 
first that a 6j symbol with fixed entries, due to the finiteness of the Racah sum rule 
(see (fT6l) in Appendix A.2), can be efficiently computed classically. On the other 
hand, the 6j is a (2d + 1) x (2d + 1) unitary matrix representing a change of basis, 

5 Recall that a quantum algorithm for solving a given computational problem is "efficient" if it 
belongs to the complexity class BQP, namely the class of problems that can be solved in polynomial 
time by a Boolean Q-circuit with a fixed bounded error in terms of the "size" of a typical input. In 
most examples the size of the input is measured by the length of the string of qubits necessary to 
encode the generic sample of the algorithmic problem, as happens with the binary representation of 
an integer number in calculations aimed to factorize it in prime factors. 
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as given explicitly in (fTTT) of Appendix A, with j'12 , j'23 representing matrix indices 
running over an interval of length 2d + 1 in integer steps. Thus the evaluation 
of the complexity class of this problem consists is asking whether, as d increases, 
the calculation of the 6j falls into the BQP class. The circuit which implements 
such task has been designed in 051 for the case of the SU(2) q 6j for each q = root 
of unity, while the analog problem involving the "classical", SU (2) 6j is still open. 

In the last few years two of the authors, in collaboration with S. Garnerone, 
have developed, on the basis of the spin network simulator setting ||29l , a new 
approach to deal with classes of algorithmic problems that classically admit only 
exponential time algorithms. The problems in questions arise in the physical con- 
text of 3D topological quantum field theories discussed in the previous section in 
the light of the fundamental result relating a topological invariant of knots, the 
Jones polynomial [36], with a quantum observable given by the vacuum expecta- 
tion value of a Wilson "loop" operator ll37l associated with closed knotted curves 
in the Witten-Reshetikhin-Turaev background model. 

Without entering into technical details, efficient (polynomial time) quantum al- 
gorithms for approximating (with an error that can be made as small as desired) 
generalizations of Jones polynomial have been found in Il35ll38ll . while the case of 
topological invariants of 3-manifolds has been addressed in j39l . The relevance in 
having solved this kind of problems stems from the fact that an approximation of 
the Jones polynomial is sufficient to simulate any polynomial quantum computa- 
tion BOH . 

Summing up, the construction of such quantum algorithms actually bears on 
the interplay of three different contexts 

1 . a topological context, where the problem is well-posed and makes it possible 
to recast the initial instance from the topological language of knot theory to 
the algebraic language of braid group theory, as reviewed in PTI ; 

2. a field theoretic context, where tools from 3D topological quantum field and 
associated 2D conformal field theory are used to provide a unitary represen- 
tation of the braid group; 

3. a quantum information context, where the basic features of quantum compu- 
tation are used to efficiently solve the original problem formulated in a field 
theoretic language. 

In the light of remark (e) at the end of section 2, further analysis of relationships be- 
tween specific 3D topological quantum field theories and (solvable) lattice models 
in 2D in the quantum-computational context would represent a major improve- 
ment not only from a theoretical viewpoint, but also in view of possible physical 
implementations. In ETl some preliminary progress has been achieved for estab- 
lishing a quantum algorithm for the evaluation of the permanent © associated with 
the partition function of the Ising model on knot-graph lattices. As shown in 11421 
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by resorting to numerical simulations, such a computational problem can be related 
to the computation of Jones invariants on suitably defined configurations, thus pro- 
viding further evidence of the "universality" of any one of the quantum algorithms 
quoted above. 

In conclusion, we hope to have been able to illustrate in sufficient details the 
role of the Wigner 6j symbol (or the q-6j) as an universal building block unify- 
ing such different fields as quantum geometry, topological quantum field theory, 
statistical lattice models and quantum computing. 

The interplay between solvability and computability within the framework of 
quantum Witten-Reshetikhin-Turaev theory and solvable lattice models deserves 
however a few more comments. Unlike perturbatively renormalizable quantum 
field theory -which represent the basic tool in the standard model in particle physics, 
where the physically measurable quantities are obtained as finite limits of infinite 
series in the physical coupling constant- quantum WRT theory is actually "solv- 
able" since functionals of type (f5]) and (O, as well as Wilson loop observables, are 
sums of a. finite number of terms for each fixed value of the deformation parameter 
q. Actually such finiteness property reflects the existence of a deeper algebraic 
symmetry stemming from braid group representations and associated Yang-Baxter 
equation, see e.g. ll37ll4Tl and references thereiij^. The issue of computability of all 
the relevant quantities of quantum WRT theory, and in particular of the Jones poly- 
nomial, is ultimately related to solvability/finiteness of the underlying theory. Thus 
the existence of "efficient" computational protocols should help in sheding light on 
the open question concerning the validation of the heuristic procedure associated 
with the path-sum quantization scheme (may be also in other contexts). Turning 
the argument upside down, the search for new efficient quantum algorithms for 
processing "invariant quantities" characterizing suitably decorated lattice, graphs, 
surfaces, etc. represents an original and possibly very fruitful approach for un- 
derstanding the underlying physical models with respect to their (yet unknown) 
integrability properties. 

Appendix A: the Wigner 6j symbol and its symmetries 

Given three angular momentum operators Ji, J2, J3 -associated with three kine- 
matically independent quantum systems- the Wigner-coupled Hilbert space of the 
composite system is an eigenstate of the total angular momentum 

Ji + J 2 + J 3 = J (V) 

and of its projection J z along the quantization axis. The degeneracy can be com- 
pletely removed by considering binary coupling schemes such as (Ji + J2) + J3 

6 This notion of solvability might be viewed as the quantum analog of the property of "complete 
integrability" in classical mechanics. Recall that integrable systems admit a sufficient number of 
conserved quantities that make it possible to solve explicitly Newton equations of motion. These 
"constants of motions" are endowed with a suitable algebraic structure under Poisson bracketing 
which is related in turn to complete integrability owing to Arnold-Liouville theorem. 
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and Ji + (J 2 + J3), and by introducing intermediate angular momentum operators 
defined by 

(Ji + J 2 ) = J12; J12 + J 3 = J (8) 



and 



(J 2 + J 3 ) = J 23 ; Ji + J23 = J, (9) 



respectively. In Dirac notation the simultaneous eigenspaces of the two complete 
sets of commuting operators are spanned by basis vectors 

\hhhih\ jm) and \j1j2hj23; jrn), (10) 

where j\ , j 2 , js denote eigenvalues of the corresponding operators, j is the eigen- 
value of J and m is the total magnetic quantum number with range —j<m< j 
in integer steps. Note that ji, j 2 , j'3 run over {0, |, 1, |, 2, . . . } (labels of SU(2) 
irreducible representations), while \ji — j 2 | < J12 < ji +32 and |j 2 — js\ < J23 < 
J2 + is (all quantum numbers are in 7i units). 

The Wigner 6j symbol expresses the transformation between the two schemes 
dHJ) and ©, namely 

\hm\23z\ jm) = V [(2 J12 + 1)(2 J23 + 1)] 1/2 {f J J f 2 ) |j'U2j3i23; j'm) 

~ [J3 J J23 J 

(ID 

apart from a phase factofl It follows that the quantum mechanical probability 

P = [(2j 12 + l)(2j 23 + l)] i jl 32 f) (12) 

U3 J J23 J 

represents the probability that a system prepared in a state of the coupling scheme 
dD, where ji,j 2 , j'3, j'12, j have definite magnitudes, will be measured to be in a 
state of the coupling scheme (O. 

The 6j symbol may be written as sums of products of four Clebsch-Gordan 
coefficients or their symmetric counterparts, the Wigner 3j symbols. The relations 
between 6j and 3j symbols are given explicitly by (see e.g. 11331 ) 

a b c\ _ \$ I ' a b c \ ( a e f ^ ( d b f\ f d e 



d e f J ^— ' /3 — 7y \a e —(pj \—S j3 tp J \S —e 7 

(13) 

where = d+e + f + 5 + e + tp. Here Latin letters stand for j'-type labels (integer 
or half-integers non-negative numbers) while Greek letters denote the associated 
magnetic quantum numbers (each varying in integer steps between — j and j, j G 
{a, b, c, d, e, /}). The sum is over all possible values of a, f3, 7, 5, e, cp with only 
three summation indices being independent. 



7 Actually this expression should contain the Racah W-coefficient W(jij2j3j; J12J23) which 
differs from the 6j by the factor Recall that (2j i2 + 1) and (2j 2 3 + 1) are the 

dimensions of the representations labeled by j'12 and j'23, respectively. 
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On the basis of the above decomposition it can be shown that the 6j symbol is 
invariant under any permutation of its columns or under interchange the upper and 
lower arguments in each of any two columns. These algebraic relations involve 
3! x 4 = 24 different 6j with the same value and are referred to as classical 
symmetries as opposite to "Regge" symmetries to be discussed in A.2. 

The 6j symbol is naturally endowed with a geometric symmetry, the tetrahe- 
dral symmetry, as the reproduction in Fig. 1 suggests. Note first that each 3j (or 
Clebsch-Gordan) coefficient vanishes unless its j-type entries satisfy the triangu- 
lar condition, namely \b — c\ < a < b + c, etc.. This suggests that each of the 
four 3j's in ([T3l can be be associated with either a 3-valent vertex or a triangle. 
Accordingly, there are two graphical representation of the 6j exhibiting its sym- 
metry properties. Here we adopt the three-dimensional picture introduced in the 
seminal paper by Ponzano and Regge [7 ], rather than Yutsis' "dual" representation 
as a complete graph on four vertices [34]. Then the 6j is thought of as a real solid 
tetrahedron T with edge lengths £i = a->r^,£2 = b+^,...,£Q = f + ^mh unit^l 
and triangular faces associated with the triads (abc), (aef), (dbf), (dec). This im- 
plies in particular that the quantities q\ = a + b + c, q2 = a + e + f,qs = b + d + f, 
q4 = c + d + e (sums of the edge lengths of each face), p\ = a + b + d + e, 
P2 = a + c + d + f,p3 = b + c + e + f are all integer with ph > qk (h = 1, 2, 3, 
k = 1, 2, 3, 4). The conditions addressed so far are in general sufficient to guaran- 
tee the existence of a non-vanishing 6j symbol, but they are not enough to ensure 
the existence of a geometric tetrahedron T living in Euclidean 3-space with the 
given edges. More precisely, T exists in this sense if {and only if, see the discus- 
sion in the introduction of [7 ]) its square volume V (T) 2 = V 2 , evaluated by means 
of the Cayley-Menger determinant, is positive. 

The features of the "quantum tetrahedron" outlined above represent the foun- 
dations of a variety of results, some of which were discovered in the golden age 
of quantum mechanics and have been widely used in old and present applications 
to atomic and molecular physics. In this paper we have tried to convey at least a 
few applications of this intriguing object in modern theoretical physics, while in 
the rest of this appendix we are going to complete the mathematical background 
needed in the previous sections, focusing in particular on semiclassical analysis 
and results from special function theory. 

A.l Ponzano-Regge asymptotic formula 

The Ponzano-Regge asymptotic formula for the 6j symbol reads [7J 

Sf^H^n)) <14) 

where the limit is taken for all entries 3> 1 (recall that h = 1) and £ r = j r + 1/2 

8 The —shift is shown to be crucial in the analysis developed in |7|: for high quantum numbers 
the length [j(j + l)] 1/2 of an angular momentum vector is closer to j + i in the semiclassical limit. 



a b d 
c f e 
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with {j r } = {a, 6, c, d, e, /}. V is the Euclidean volume of the tetrahedron T and 
9 r is the angle between the outer normals to the faces which share the edge £ r . 
From a quantum mechanical viewpoint, the above probability amplitude has the 
form of a semiclassical (wave) function since the factor 1 / \/24irV is slowly vary- 
ing with respect to the spin variables while the exponential is a rapidly oscillat- 
ing dynamical phase. Such kind of asymptotic behavior complies with Wigner's 
semiclassical estimate for the probability, namely { a c h j e } ~ 1/12-7T V , to be 
compared with the quantum probability given in (PT2l) . Moreover, according to 
Feynman path sum interpretation of quantum mechanics R3l . the argument of the 
exponential in (fl4l) must represent a classical action, and indeed it can be read 
as ^ pq for pairs (p,q) of canonical variables (angular momenta and conjugate 
angles). Such an interpretation has been improved recently by resorting to multidi- 
mensional WKB theory for integrable systems and geometric quantization methods 

ma. 





■ a p 


b x .. 


■ b q 



A.2 Racah hypergeometric polynomial 

The generalized hypergeometric series, denoted by p F q , is defined on p real or 
complex numerator parameters a\, 02, ■ ■ ■ , a p , q real or complex denominator pa- 
rameters b\ , 62 , ■ ■ ■ , b q and a single variable z by 

j Qp . \ _ (ai)n---(ap)n z n 5 

where (a) n = a(a+l)(a+2) ■ ■ ■ (a+n—1) denotes arising factorial with (a)o = 1. 
If one of the numerator parameter is a negative integer, as actually happens in the 
following formula, the series terminates and the function is a polynomial in z. 
The key expression for relating the 6j symbol to hypergeometric functions is given 
by the well-known Racah sum rule (see e.g. 0, topic 11 and Il33l . Ch. 9 also 
for the original references). The final form of the so-called Racah polynomial is 
written in terms of the 4F3 hypergeometric function evaluated at z = 1 according 
to 

" b f d \ = A(abe)A(cde)A(acf)A(bdf) + 1)! 



, ai—Pi ct2—f3\ Q3— /3i 04—^1 



where 



(/3 2 - hWz - P1W1 - ai)!(/3i - a 2 )!(/3i - /3 3 )!(A - a 4/ 



f3\ = min(a + b + c + d; a + d + e + / ';b + c + e + /) 



(16) 



and the parameters 02 , P3 are identified in either way with the pair remaining in the 
3-tuple {a + b + c + d;a + d + e + f]b + c + e + f) after deleting j3\. The four a's 
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may be identified with any permutation of (a + b+e; c+d+e; a + c+f; b+d+f). 
Finally, the A-factors in front of 4F3 are defined, for any triad (abc) as 



\ / 1 \ \{a + b-c)\{a-b + c)\{-a + b + c)\V /2 
A (abc) = -. : -r- 

Such a seemly complicated notation is indeed the most convenient for the purpose 
of listing further interesting properties of the Wigner 6j symbol. 

• The Racah polynomial is placed at the top of the Askey hierarchy including 
all of hypergeometric orthogonal polynomials of one (discrete or continuous) 
variable [45]. Most commonly encountered families of special functions 
in quantum mechnics are obtained from the Racah polynomial by applying 
suitable limiting procedures, as recently reviewed in li46l . Such an unified 
scheme provides in a straightforward way the algebraic defining relations of 
the Wigner 6j symbol viewed as an orthogonal polynomial of one discrete 
variable, cfr. (fT6l ). By resorting to standard notation from the quantum theory 
of angular momentum, the defining relations are: 

the Biedenharn-Elliott identity (R = a + b + c + d + e + f + p->rq + r): 



p q r I J p q r 
e a d( \f b c 



(17) 



the orthogonality relation (5 is the Kronecker delta) 



• Given the relation (TT6l) . the unexpected new symmetry of the 6j symbol dis- 
covered in 1958 by Regge [47] (see also (T][33l) is recognized as a "trivial" 
set of permutations on the parameters a, /3 that leaves 4F3 invariant. Com- 
bining the Regge symmetry and the "classical" ones, one get a total number 
of 144 algebraic symmetries for the 6j. Note however that implications of 
Regge symmetry on the geometry of the quantum tetrahedron, taken into 
account in 11481 . certainly deserve further investigations also in view of the 
relevance of this topic in completely different contexts, cfr. for instance [49]. 

• The Askey hierarchy of orthogonal polynomials can be extended to a q- 
hierarchy [ 431 . on the top of which the q-4-F3 polynomial stands. 

It is worth noting that the deformation parameter q was originally assumed 
by physicists to be a real number related to Planck constant h by q = e h , and 
therefore it is commonly referred to as a 'quantum' deformation, while the 
'classical', undeformed Lie group symmetry is recovered at the particular 
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value q = 1. However, when dealing with quantum invariants of knots and 
3-manifolds formulated in the framework of "unitary" quantum field theory, 
as done in section 2 and 3, q is taken to be a complex root of unity, the case 
q = 1 being considered as the "trivial" one. We refer to |[50ll5Tl for accounts 
on the theory of q-special functions and q-tensor algebras. 
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